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Abstract 

The elliptic-flow ratio of neutrons with respect to protons or light complex par- 
ticles in reactions of neutron-rich systems at relativistic energies is proposed as an 
observable sensitive to the strength of the symmetry term in the equation of state at 
supra-normal densities. The results obtained from the existing FOPI/LAND data 
for 197 Au + 197 Au collisions at 400 MeV/nucleon in comparison with the UrQMD 
model favor a moderately soft symmetry term with a density dependence of the 
potential term proportional to (p/po)' y with 7 = 0.9 ± 0.4. 
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Recent advances in the observation and modelling of astrophysical phenom- 
ena have renewed the interest in the nuclear equation of state (EOS) and, in 
particular, in its dependence on density and on asymmetry, i.e., on the rel- 
ative neutron-to-proton abundance |T1[2l[3] . Supernova simulations or neutron 
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star models require inputs for the nuclear EOS at extreme values of these pa- 
rameters [US]. The dependence on asymmetry is given by the symmetry term 
whose evolution with density is not only of interest for astrophysics but also 
for nuclear physics. The thickness of the neutron skin of heavy nuclei, e.g., 
reflects the differential pressure exerted on neutrons in the core [6], and the 
strength of the three-body force, an important ingredient in nuclear structure 
calculations [TJ, represents one of the major uncertainties in modelling the 
EOS at high density [PIS] . 

Considerable efforts are, therefore, underway in using heavy-ion reactions for 
extracting experimental information on the symmetry energy which is the dif- 
ference between the energy of neutron matter and of symmetric matter. While 
fairly consistent constraints for the symmetry energy near normal nuclear mat- 
ter density have been deduced from recent data [31I9"1I10|I11] , much more work is 
still needed to probe its high-density behavior. The predictions of microscopic 
models diverge widely there [3~fT2] . and the interpretation of obvious observ- 
ables turns out to be more difficult than perhaps anticipated. The ratios of 
K + /K° production were found to be less sensitive to the symmetry energy 
when actual collisions were modelled rather than matter under equilibrium 
conditions |l3j . Variations of the tc~ /tt + ratio of up to 20% for soft versus stiff 
parameterizations are expected but the competition of mean-field and collision 
effects appears to be very delicate. Different predictions leading to practically 
opposite conclusions were obtained from the comparison of the FOPI data for 
197 Au + 197 Au collisions [TJ] with different transport models [T5"fl6|ll7] . 

It seems essential, in this situation, to enlarge the experimental basis of suit- 
able probes and to extend it from isotopic yield ratios to the isospin depen- 
dence of dynamical observables. Densities of up to 2 or 3 times the saturation 
density may be reached within a short time scale (~ 20 fm/c) in the central 
zone of heavy-ion collisions in the present range of relativistic energies of up 
to ~ 1 GeV/nucleon [18]. The resulting pressure produces a collective out- 
ward motion of the compressed material whose strength will be influenced by 
the symmetry energy in asymmetric systems [19]. Flow observables have been 
proposed by several groups as probes for the equation of state at high den- 
sity [2"0"f2"Tj . among them the so-called differential neutron-proton flow which 
is the difference of the parameters describing the collective motion of free neu- 
trons and protons weighted by their numbers [20]. As pointed out by Yong 
et al. [22], this observable minimizes the influence of the isoscalar part in the 
EOS while maximizing that of the symmetry term. Its proportionality to the 
particle multiplicities, however, makes its determination very dependent on 
the experimental efficiencies of particle detection and identification and on 
the precise theoretical distinction between free and bound nucleons. In this 
work, the ratios of neutron versus proton or light-charged-particle flows will 
be considered and their sensitivities to the strength of the symmetry term will 
be explored. The problem of systematic errors is thereby reduced. 
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Fig. 1. Scatter plot of 10000 neutron events with energy > 40 MeV from a run 
without shadow bars in the plane of transverse momentum pt vs. rapidity y. The 
quality criteria are the same as used in the flow analysis. Their application produces 
the local inefficiency in the forward detector unit caused by slabs with reduced 
performance. The dashed lines represent the acceptance cuts pt > 0.3 GeV/c and 
0.25 < y/y p < 0.75 applied later in the analysis. 

In a series of experiments at the GSI laboratory combining the LAND and 
FOPI (Phase 1) detectors, both neutron and hydrogen collective flow observ- 
ables from 197 Au + 197 Au collisions at 400, 600 and 800 MeV/nucleon have 
been measured [23f24~] . These data sets are presently being reanalyzed in or- 
der to determine optimum conditions for a dedicated new experiment, but 
also with the aim to produce constraints for the symmetry energy by compar- 
ing with predictions of state-of-the-art transport models. Here, we report the 
results obtained with the 400- MeV/nucleon data set and the Ultrarelativistic- 
Quantum- Molecular-Dynamics (UrQMD) model which has recently been adapted 
to heavy ion reactions at intermediate energies [25] . 

The Large- Area- Neutron-Detector LAND is a 2x2x1 m 3 calorimeter consisting 
of in total 200 slabs of interleaved iron and plastic strips viewed by photomul- 
tiplier tubes at both ends [26J. For the flow experiment, LAND was divided 
into two parts of 50-cm depth each, positioned to cover laboratory angles of 
37°-53° and 61°-85°. The corresponding acceptance in the plane of transverse 
momentum p t vs. rapidity y is shown in Fig. [TJ Veto walls of 5-mm-thick 
plastic scintillators in front of the two detector units permitted the distinc- 
tion of neutral and charged particles and to resolve atomic numbers Z < 2. 
A shower-recognition algorithm is used to identify individual particles within 
the registered hit distributions. Particle momenta were determined from the 
measured time-of-flight (At « 550 ps, FWHM) over flight paths of 5 m and 
7 m to the more backward and more forward detector units, respectively. The 
amount of light collected from identified showers was used to generate mass 
spectra of the hydrogen isotopes as described in Ref. [24J. Individual isotopic 



resolution is not achieved but the strong proton group is clearly identified. The 
selected proton samples, containing about 45% of the total hydrogen yields, 
are expected to have isotopic purity of better than 95%. 

During its Phase 1, the 4-ir detector FOPI consisted of a forward plastic array, 
built from more than 700 individual plastic-scintillator strips and covering 
the forward hemisphere at laboratory angles from 1° to 30° [27J. LAND was 
used to detect and identify neutrons and light charged particles with nearly 
identical methods while the FOPI Forward Wall was used to determine the 
modulus and orientation of the impact parameter from the multiplicity M c 
and distribution of the detected charged particles [28]. The background of 
scattered neutrons was measured in separate runs with shadow bars made 
from iron and positioned to block the direct flight paths from the target to 
one or the other of the two LAND units. Additional details can be found in 
Refs. [25123] . 

The usefulness of the studied observables has been evaluated with the UrQMD 
transport model [25]. This model, originally developed to study particle pro- 
duction at high energy [29], has been adapted to intermediate-energy heavy- ion 
collisions by introducing a nuclear mean field corresponding to a soft EOS with 
momentum dependent forces and with different options for the dependence on 
asymmetry. Two of these choices are used here, expressed as a power-law de- 
pendence of the potential part of the symmetry energy on the nuclear density 
p according to 

E sym = + JSgS, = 22 MeV • (p/ Po ) 7 + 12 MeV • (p/p ) 2/3 (1) 

with 7 = 0.5 and 7 = 1.5 corresponding to a soft and a stiff density dependence 
(po is the normal nuclear density of ~ 0.16 fm -3 ). The kinetic part remains 
unchanged. The UrQMD model has been successfully used for studies of a wide 
range of heavy-ion-collision problems at intermediate and relativistic energies 
(Ref. [30] and references given therein). 

As mentioned above, a realistic description of the clustering processes during 
the temporal evolution of the reaction is crucial for predicting dynamical prop- 
erties of free neutrons, protons and light charged particles. In the UrQMD, the 
clustering algorithm is based on the proximity of nucleons in phase space with 
two parameters for the relative nucleon coordinates and momenta. The present 
analysis was performed with the distributions obtained after a reaction time 
of 150 fm/c and with the proximity limits Ar = 3.0 fm and Ap = 275 MeV/c 
which are typical for QMD models [3T] . These parameters had been slightly 
adjusted from their default values in order to obtain a best possible descrip- 
tion of measured fragment spectra as a function of Z. The result for central 
collisions of 197 Au + 197 Au at 400 MeV/nucleon is shown in Fig. [2] in com- 
parison with the data of Reisdorf et al. [25]. With a normalization at Z = 1, 
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Fig. 2. Fragment yields, integrated over the 4-7r solid angle, in central (equivalent 
to b < 2.0 fm) collisions of 197 Au + 197 Au at 400 MeV/nucleon as a function of Z 
(dots, from Ref. [28J) in comparison with UrQMD predictions normalized at Z = 1 
(histogram). 

the overall dependence on Z is rather well reproduced but the yields of Z = 2 
particles are underpredicted by about a factor of 3. The strong binding of 
a particles is beyond the phase-space criterion used in the model. However, 
also the 4-7r integrated yields of deuterons and tritons in central collisions are 
underestimated by similar factors of 2 to 3. 

Flow observables were obtained from the calculated event samples by fitting 
the azimuthal particle distributions with the usual Fourier expansion 

/(A0) oc 1 + 2 • V! • cos(A0) + 2 ■ v 2 ■ cos(2 ■ A<p) (2) 

with A(j) representing the azimuthal angle of the momentum vector of the 
emitted particle with respect to the reaction plane [32]. The predictions ob- 
tained for the elliptic flow of neutrons, protons, and hydrogen yields for 197 Au 
+ 197 Au at 400 MeV/ nucleon and for the two choices of the density dependence 
of the symmetry energy, labeled asy-stiff (7 = 1.5) and asy-soft (7 = 0.5), are 
shown in Fig. [3l The dominant difference is the significantly larger neutron 
squeeze-out in the asy-stiff case (upper panel) compared to the asy-soft case 
(lower panel). The proton and hydrogen flows respond only weakly, and in 
opposite direction, to the variation of 7 within the chosen interval. 

Collective flow has recently been measured rather systematically |33j . The re- 
liability of the methods has been demonstrated by the good agreement of data 
sets from different experiments in the overlap regions of the studied intervals 
in collision energy [32] . An excitation function of flow for 197 Au + 197 Au col- 
lisions has been produced for incident energies from below 100 MeV/nucleon 
up to the ultrarelativistic regime. Within this wide range, the negative ellip- 
tic flow of light particles, i.e. their squeeze-out perpendicular to the reaction 
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Fig. 3. Elliptic flow parameter V2 for mid-peripheral (5.5 < b < 7.5 fm) 197 Au + 
197 Au collisions at 400 MeV/nucleon as calculated with the UrQMD model for neu- 
trons (dots), protons (circles), and all hydrogen isotopes (Z = 1, open triangles), 
integrated over transverse momentum pt, as a function of the laboratory rapidity 
yiab- The predictions obtained with a stiff and a soft density dependence of the sym- 
metry term are given in the upper and lower panels, respectively. The experimental 
result from Ref. [33] for Z = 1 particles at mid-rapidity is represented by the filled 
triangle (the horizontal bar represents the experimental rapidity interval). 



plane, has its maximum at 400 MeV/nucleon. The value v?. = —0.082 mea- 
sured at this energy for Z = 1 particles at mid-rapidity [33] is reasonably well 
reproduced by the UrQMD model, even though underestimated by 8% and 
15% with the two parameterizations (Fig. [3]). 

The predictions for elliptic particle flows depend strongly on the type of 
transport theory used and on its parameters. Values lower than found ex- 
perimentally have been reported also for other models as, e.g., discussed in 
Refs. [3Tf34"] . In a recent paper, Li et al. investigate the elastic in- medium 
nucleon-nucleon cross section and show the effects of its momentum depen- 
dence on various observables [35]. The standard parameterization, also used 
here up to now, is labeled FP1 in that paper. One of the alternatives, FP2, 
exhibits a higher correction factor than FP1 for relative nucleon-nucleon mo- 
menta around 0.2 < p NN < 0.5 GeV/c at high density. With this parameteri- 
zation, the elliptic flow is larger, now exceeding the measured Z = 1 value by 
16% and 19% for the stiff and soft parameterizations of the symmetry term, 
respectively. 
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Fig. 4. Measured flow parameters v\ (top) and V2 (bottom) for mid-peripheral 
(5.5 < b < 7.5 fm) 197 Au + 197 Au collisions at 400 MeV/nucleon for neutrons 
(dots), protons (circles), and hydrogen isotopes (Z = 1, open triangles) integrated 
within 0.3 < pt/A < 1.3 GeV/c/nucleon as a function of the normalized rapidity 
yiab/lJp- The UrQMD predictions for neutrons, obtained with a stiff (7 = 1.5, full 
lines) and a soft (7 = 0.5, dashed) density dependence of the symmetry term have 
been filtered to correspond to the geometrical acceptance of the experiment. The 
experimental data have been corrected for the dispersion of the reaction plane (see 
text). 

The results obtained from the FOPI/LAND data set for 400 MeV/nucleon 
are shown in Fig. H as a function of the rapidity y, normalized with re- 
spect to the projectile rapidity y p = 0.896, and for the PM3 event class 
(27 < M c < 39) which approximately corresponds to the interval of impact 
parameters 5.5 fm < b < 7.5 fm [28]. The geometrical acceptance of the LAND 
detector in this experiment produces a dependence of the accepted range in 
transverse momentum p t on rapidity y, pt increases with y (Fig. [T|), and is re- 
sponsible for the observed asymmetry of v 1 and V2 with respect to mid-rapidity 
y/Up — 0-5. Gates were set for transverse momenta 0.3 < Pt/A < 1.3 GeV/c, 
mainly to ensure equal lower thresholds for neutrons and charged particles 
detected with LAND. For the comparison, the UrQMD outputs have been 
filtered with the geometrical acceptance of the LAND detector and subjected 
to the same transverse-momentum gates as the experimental data. The ex- 
perimental Fourier coefficients, on the other hand, have been corrected for the 
uncertainty associated with the reconstruction of the reaction plane from the 
azimuthal distribution of particles measured with the FOPI Forward- Wall de- 
tector. Because of the significant multiplicities M > 26 of particles observed 
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in the considered range of impact parameters b < 7.5 fm, this correction is 
not crucial [32], and global correction factors 1.05 for V\ and 1.15 for t> 2 have 
been adopted. 

Only the predictions for neutrons, obtained with the standard FP1 parame- 
terization of the elastic nucleon-nucleon cross section [35] , are shown in Fig. HI 
for clarity as well as because of the much smaller dependence on 7 of the 
results for hydrogen isotopes (cf. Fig. [3]). A nearly negligible sensitivity to the 
stiffness of the symmetry energy is exhibited by the directed flow, according to 
the UrQMD model (Fig. HI top panel). The lines representing the predictions 
for the soft and stiff parameterizations of the symmetry energy fall practically 
on top of each other. The results for protons or hydrogen isotopes are nearly 
identical. All predictions compare well, however, with the experimental results 
for the multiplicity bin PM3 corresponding to this range of mid-peripheral im- 
pact parameters. The range of rapidities y > 0.8 at which the data deviate 
from the expected linearity coincides with the shift of the LAND acceptance 
to transverse momenta pt > 0.75 GeV/c/nucleon at which the yields start to 
drop rapidly (Fig. [TJ). 

Also the second Fourier coefficient v 2 describing elliptic flow and its depen- 
dence on rapidity are fairly well described by the UrQMD calculations (Fig. HJ 
bottom panel). Here the statistical errors of the experimental data are larger 
but the slight underprediction expected from Fig. [3] is clearly recognized. Con- 
trary to the directed flow, the v 2 values calculated with the asy-stiff and asy- 
soft parameterizations for neutrons are significantly different from each other. 
This is not the case for the hydrogen isotopes whose squeeze-out parameters 
are close to those of neutrons for the asy-stiff case or, for free protons alone, 
similar to the asy-soft neutron predictions. Even though the experimental er- 
rors are considerable, the sensitivity of the neutron squeeze-out seems large 
enough to permit the extraction of a useful result from the present data set. 

The dependence of the elliptic flow parameter v 2 on the transverse momentum 
per nucleon, Pt/A, is shown in Fig. [5j upper panel, for the combined data set 
of central and mid-peripheral collisions (PM3-PM5, M c > 27, cf. Ref. [28J) 
in order to exploit the full statistics collected in the experiment. The extreme 
rapidities are avoided by selecting the interval 0.25 < y/y p < 0.75. The in- 
crease of v 2 in absolute magnitude is nearly linear and the measured values 
for neutrons and hydrogen isotopes follow approximately the UrQMD predic- 
tions for b < 7.5 fm. For hydrogens (not shown), the results for the stiff and 
soft density dependences are similar and close to the asy-stiff predictions for 
neutrons. For neutrons, however, the predictions are different by, on average, 
20% for the main part of the pt interval (dashed lines) but seem to converge 
at high p t at which the yields become small and the statistical errors large. 

For the quantitative evaluation, the ratio of the flow parameters of neutrons 
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Fig. 5. Differential elliptic flow parameters t>2 for neutrons (dots) and hydrogen 
isotopes (open triangles, top panel) and their ratio (lower panels) for moderately 
central (b < 7.5 fm) collisions of 197 Au + 197 Au at 400 MeV/nucleon, integrated 
within the rapidity interval 0.25 < y/y p < 0.75, as a function of the transverse 
momentum per nucleon pt/A. The symbols represent the experimental data. The 
UrQMD predictions for 7 = 1.5 (a-stiff) and 7 = 0.5 (a-soft) obtained with the 
FP1 parameterization for neutrons (top panel) and for the ratio (middle panel), 
and with the FP2 parameterization for the ratio (bottom panel) are given by the 
dashed lines. Note the suppressed zero of the abscissa. 

versus hydrogen isotopes is proposed to be used. Systematic effects influencing 
the collective flows of neutrons and charged particles in similar ways will be 
minimized in this way, on the experimental as well as on the theoretical side. 
Among them are, in particular, the still existing uncertainties of isoscalar 
type in the EOS but also, e.g., the dispersion of the experimentally deter- 
mined orientation of the reaction plane, possible small azimuthal anisotropies 
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in the particle detection with the Forward Wall, and the matching of the 
impact-parameter intervals used in the calculations with the experimental 
event groups selected according to multiplicity. 

The results for the ratio v^/v^, i.e. with respect to the integrated hydrogen 
yield, are shown in Fig. [5] (lower panels). They exhibit once more the sensitivity 
of the elliptic flow to the stiffness of the symmetry energy predicted by the 
UrQMD. It is slightly smaller with the FP2 parameterization of the in-medium 
nucleon-nucleon cross section (bottom panel). The experimental ratios, even 
though associated with large errors, are found to scatter within the intervals 
given by the calculations for 7 = 0.5 and 1.5 in either case, indicating that 
the effects of the cross section parameterization largely cancel in the ratios. 
Linear interpolations between the predictions, averaged over 0.3 < pt/A < 1.0 
GeV/c with weights (A7j)~ 2 derived from the experimental errors A^i for 
each p t bin, yield 7 = 1.01 ± 0.21 and 7 = 0.98 ± 0.35 for the FP1 and FP2 
parameterizations, respectively. 

The same analysis was also performed for the squeeze-out ratios v^jv^ of 
neutrons with respect to free protons. The resulting power-law exponents are 
7 = 0.99 ± 0.28 and 7 = 0.85 ± 0.47 for the FP1 and FP2 parameterizations, 
respectively. They are nearly identical to those obtained from the neutron-to- 
hydrogen squeeze-out ratios, even though associated with larger statistical un- 
certainties. If the comparisons of the neutron-to-hydrogen ratios v^/v^ are re- 
stricted to the mid-peripheral interval of impact-parameters 5.5 < b < 7.5 fm, 
smaller values 7 = 0.58 ± 0.27 and 7 = 0.35 ± 0.44 are obtained with the 
FP1 and FP2 parameterizations, respectively. This apparent dependence on 
impact parameter existed already in preliminary results reported earlier [36]. 

All results are, within errors, compatible with each other. In particular, the fact 
that the yields of deuterons and tritons with respect to protons are underesti- 
mated in the model calculations does not seem to have visible consequences. 
The ratios with respect to all hydrogen isotopes and to protons alone yield 
the same exponent 7. However, because of the limited statistical accuracy, the 
existence of systematic uncertainties cannot be ruled out and may, in fact, be 
partly responsible for the apparent dependence on impact parameter. We will, 
therefore, assume that the difference between the weighted means of the four 
results for b < 7.5 fm and the full statistics, 7 = 0.98, and of the two values 
for V2/V2 and the mid-peripheral range 5.5 < b < 7.5 fm, 7 = 0.52, can serve 
as an estimate, A7| syst . = ±0.23, for the magnitude of possible systematic 
uncertainties. 

The statistical errors of the elliptic flow coefficients translate into statisti- 
cal errors of 7 that depend on the sensitivity displayed by the predictions 
(Fig. [5l lower panels). Since the magnitude of elliptic flow is either slightly 
underpredicted or overpredicted with the two parameterizations, we interpo- 
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late between the two cases and obtain A7| stat . = ±0.27 as the most probable 
statistical error. Combining it with the systematic uncertainty and considering 
the spread and trends of the individual results, we adopt the value 7 = 0.9±0.4 
as the result of this analysis for the exponent describing the density depen- 
dence of the potential term. Together with the kinetic term proportional to 
(p/po) 2 / 3 , the squeeze-out data thus indicate a moderately soft behavior of 
the symmetry energy. It is obvious that data with higher statistical precision 
are highly desirable. Being able to follow the evolutions with impact parame- 
ter, transverse momentum (cf. Fig. [5]) and particle type more accurately will 
provide supplementary constraints for the comparison with theory and will 
permit a clearer identification of systematic dependences related to it. 

Other systematic errors were found to be comparatively minor. In the ex- 
perimental data analysis, the explored variations of the conditions and gates 
within reasonable limits have only small effects. The uncertainty associated 
with the procedures used to measure and subtract the background of scattered 
neutrons detected with LAND |23|24] has been estimated by subtracting only 
60% instead of 100% of the measured background contributions from the total 
neutron yields. It reduces the resulting power-law exponent 7 by IA7I « 0.2 
which may be considered an upper limit for systematic experimental effects. 

Additional systematic uncertainties on the theoretical side have also been eval- 
uated. In particular, a test has been made whether supra-saturation densities 
are actually probed by the squeeze-out ratios. The coefficient of the potential 
term in the symmetry energy has been reduced from its default value 22 MeV 
to 18 MeV (cf. Eq.flTJ)). A softer (stiffer) density dependence is expected in 
this case if the symmetry term were mainly tested at densities below (above) 
saturation. The result obtained for v^jv\ with FP1, 7 = 0.93 ± 0.28, is close 
to the corresponding 7 = 1.01 ± 0.21 obtained with 22 MeV, thus confirming 
that the dynamics at supra-saturation density are important for the compar- 
ison. The parameter L, proportional to the slope of the symmetry term at 
saturation [3] changes by less than 20 MeV when this modification is made. 
The adopted result 7 = 0.9 ± 0.4 corresponds to L = 83 ± 26 MeV in the 
standard parameterization with 22 MeV. 

In summary, the ratio of the elliptic-flow parameters of neutrons with re- 
spect to hydrogen isotopes measured for 197 Au + 197 Au collisions at 400 
MeV/nucleon in the FOPI/LAND experiment has been used to determine 
the strength of the symmetry term in the equation of state by comparing it 
to UrQMD model predictions. A value of 7 = 0.9 ± 0.4 has been obtained 
for the power-law coefficient describing the density dependence of the poten- 
tial term. Individual systematic errors were not found to exceed a margin of 
IA7I rs 0.2. Together with the kinetic term proportional to (p/po) 2 ^ 3 , this 
indicates a moderately soft symmetry term. Within errors, it is consistent 
with the density dependence deduced from fragmentation experiments prob- 
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ing nuclear matter near or below saturation [9|10|lll] but inconsistent with the 
super-soft or extremely stiff behavior obtained from the comparisons of the 
tt~/tt + yield ratios measured for the same reaction with different transport 
models jTrJyTT] . Further experimental and theoretical work will be needed to 
resolve these apparent ambiguities. 
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